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ABSTRACT 


In this report we present the results of our research activities during the past year 
(FY 1 972-1973) on holographic grating research. A large portion of this work was performed 
using rigorous vector diffraction theory; therefore, we have included the necessary theory in 
this report. The diffraction efficiency studies begun in the first year (FY 1971-1972) were 
continued using programs based on a rigorous theory. The simultaneous occurrence of high 
diffraction efficiencies and the phenomenon of double Wood’s anomalies is demonstrated 
along with a graphic method for determining the necessary grating parameters. Also, an 
analytical solution for a grating profile that is perfectly blazed is obtained. The performance 
of the perfectly blazed grating profile is shown to be significantly better than grating pro- 
files previously studied. Finally, we describe a proposed method for the analysis of coarse 
cchellc gratings using rigorous vector diffraction that is currently being developed. 
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1. INTRODUCTION 

During flic past two years our research activities have been directed at the analysis, de- 
sign, and evaluation of diffraction gratings produced by holographic techniques. During the 
lirsl year of this contract (1971-1972) we applied approximate equations, derived from vec- 
tor diffraction theory, to the study of diffraction gratings with different groove profiles. In 
particular, we studied echelette profiles that have rounded edges. This profile is very similar 
to those produced using holographic methods. For these studies computer programs, based 
on the approximate theory, were developed and tested to determine their regions of valid- 
ity. The programs were then used to study the efficiency characteristics to be expected from 
holographic gratings. We were able to show that holographic gratings can be designed so that 
nearly all of the incident radiant energy is contained in a preselected diffraction order for a 
particular wavelength. It was also shown that these high blaze efficiencies occur in the pres- 
ence of a double Wood’s anomaly. 

During the following year of this contract (1972-1973) we continued the above studies. 
The earlier studies had shown that the computer programs based on the approximate theory 
worked very well for gratings with shallow grooves but that the approximations were not 
valid for gratings with deep grooves. To avoid a duplication of effort, we obtained a set of 
computer programs based on rigorous vector diffraction theory and valid for deep echelette 
gratings. These programs were made available by Professor A. R. Nereuther at the University 
of California, Berkeley. The programs are based on the works of Dr. K. A. E. H. Zaki 1 and 
Dr. H. A. Kalhor. 2 Throughout this report the results obtained from the use of these pro- 
grams will be clearly distinguished from those obtained using our own programs. 

The Berkeley programs were first modified so that they were easier to use. Also, an 
attempt was made to modify the programs so that they could analyze gratings with an arbi- 
trary profile. However, this effort was not successful as the results obtained from the pro- 
grams were no longer always reliable. Thus, the latter modifications were removed from the 
programs. 

The study of the design of blazed diffraction gratings was extended to include an analyti- 
cal solution to the problem. It was shown that it was possible to design a grating profile so 
that 1 00% of the incident energy was contained in a preselected diffraction order. In addi- 
tion, the grating was now well blazed over a broad wavelength interval as opposed to a single 
wavelength as was obtained in the previous study mentioned earlier. 

It became apparent, through close contact with our contract monitor and his associates, 
that the needs of NASA-GSFC far exceeded the capabilities of the programs based on the 
approximate theory as well as the Berkeley programs in their present form. This was pri- 
marily due to the very deep gratings and the large number of real diffraction orders resulting 
from coarse echelle gratings that were used in the vacuum ultraviolet. In order to analyze 
these gratings we directed all of our recent efforts toward the formulation of an efficient 



method by which the equations derived from rigorous vector diffraction theory could be 
applied to the analysis of coarse echelle gratings as required by NASA. The salient features 
of the rigorous theory are discussed in the next section. In Section 5 we outline a method 
that can be used for the study of coarse echelle gratings. A more complete mathematical 
derivation of the equations required for the study is presented in Appendix B, C, and D. 
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2. VECTOR DIFFRACTION FROM DEEP GROOVE 
DIFFRACTION GRATINGS 

The grooves in a diffraction grating arc said to be deep when the groove depth approaches 
or is greater than the wavelength of the incident radiant energy. To analyze such a grating a 
rigorous theory that allows the amplitude coefficients of the diffracted fields to be a func- 
tion of position within the grating profile must be used. It is apparent that the assumption 
that the amplitude coefficients are constant within the grating groove will become less valid 
as the grating depth increases. In fact, this is exactly why the approximate theory fails when 
used to analyze deep gratings. Below we describe a rigorous theory that can be applied to 
the analysis of deep echelle gratings. 

Consider a perfectly conducting, planar reflection grating that is irradiated by a mono- 
chromatic plane wave. We wish to determine the amount of energy appearing in each of the 
diffracted orders. This is done by using Maxwell’s electromagnetic field equations and apply- 
ing the appropriate boundary conditions. 

The total electromagnetic field is defined as the incident field plus the diffracted fields 

Total field = (Incident field) + (Diffracted fields). 

This applies to both the electric field E as well as the magnetic field H. We assume a right- 
handed coordinate system such that the y axis is perpendicular to the plane of the grating, 
the x axis is perpendicular to the grating grooves, and the 2 axis is parallel with them. The 
grating profile is given by 

y = fix). (1) 

An electromagnetic wave with any arbitrary state of polarization can be decomposed into 
two waves linearly polarized at right angles to each other. Therefore, we consider two cases: 
In the first the electric field vector is parallel to the grating grooves and is called the E- 
parallel state of polarization. In the second the magnetic field vector is parallel to the grating 
grooves and is called the //-parallel state of polarization. This decomposition reduces the 
vector equations to scalar equations. On the surface of the grating the boundary conditions 
to be obeyed by each state of polarization are given by 

(a) /f-parallel 

£ T oT(*-T)L n s ^ 

b=/U) 

That is, the tangential component of the total electric field on the surface of the grating 
must be zero. 
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(b) //-parallel 


AZ/tot^’^) 

dn 


= 0 . 

v = /W 


( 3 ) 


which says that the normal derivative of the total magnetic field on the surface of the grat- 
ing must be zero. 

With the use of Maxwell’s equations and the application of the above boundary condi- 
tions it is possible to calculate the desired diffraction efficiencies. Below we have briefly 
described the procedure used to perform this calculation. The procedure we describe is 
exactly the same as that used by Zaki and Kalhor. 

The first step is to calculate the field due to an infinitely long alternating line current by 
using Maxwell’s equations. This information is then used to determine the fields due to an 
infinite array of equally spaced line currents. The result of this calculation is then used to 
formulate a suitable Green’s function that describes the problem. In formulating the Green’s 
function the currents are assumed to be excited by a plane wave incident at an angle 6 t . The 
Green’s function so obtained is a function of x,y and x ' ,y' separately. That is, the Green’s 
function describes the fields at the point x,y due to an infinite line current at the point 
x'V. Due to the choice of the phasing of the individual line currents, the Green’s function 
exhibits the following symmetry 

G(x+ma,y;x'+pa,y') = exp(—ikma sin0 ( ) G{x,y,x ,y‘) (4) 

where a is the grating constant and m and p are positive or negative integers. We note that 
for a normally incident plane wave ((?,• = 0) the Green’s function is strictly periodic in both 
variables x and 

Having obtained the Green’s function we then use Green’s identity to obtain the dif- 
fracted fields. The surface integral that appears in Green’s identity is evaluated by consider- 
ing the volume enclosed by the surfaces shown in Fig. 1 . The surface -S', is along the grating 
surface. The surfaces S 5 and S 6 (not shown in the figure) are parallel to the plane of the 
paper and are the front and rear surfaces of the closed volume. It can be shown that the 
contributions due to S s and S 6 are equal and opposite as are the contributions due to S 2 
and S 3 . The contribution due to the integration along S A can easily be determined. There- 
fore, we find that the field at any point within the volume that is enclosed by the different 
surfaces is given by the surface integral along the grating profile. For the ZT -parallel polariza- 
tion we obtain 


£c/{-v ,y ) - exp{+ikx sinfl,) f G(x,y;x ,y ) — ~ — dl 

J c on 


( 5 ) 




where E d (x',y') is the diffracted field, and if EJx.y) is the incident field then the total 
electric field £(x,jO is given by 

E(x,y) = Ej(x,y) + E d (x,y). (6) 

The integration is along the grating profile with the normal directed into the volume as 
shown in Fig. 2. 



Fig. 2. The surface normal is directed into the vol- 
ume enclosed by the surface used to evaluate the 
surface integral. 
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For the //-parallel polarization we obtain the result 


H d {x’y ) = cxp(+ikx sin0 ; ) £ £//(*. j) — A - - j dl ( 7 ) 

As was the case with Eq. (5), H d is the diffracted field, and the total field H is given by the 
sum of the incident field plus the diffracted fields 

IHx.y) = H,(x,y) + H d (x,y). {8) 

In Eqs. (5) and (7) the elementary length {dl) for the contour integral along the grating 
profile may be converted to an integral along x by observing that 



The limits of integration for the variable X are from jc = 0 to x = a. 

It is also useful to note that the normal derivative is given by 

A = r 1+ (my>_YY' A r (dm\ a , a i 

L \ dx ) J [ \ 3jc / dx by J ' (10) 

We may now write Eqs. (5) and (7) in a more convenient form by adding the respective 
incident fields to both sides and invoking the boundary conditions on the surface of the 
grating. For the ^-parallel case we find that 


Ei(x'.J(x')) = + exp{ikx' sin0 r ) j G(x,f{x);x\f{x))^{ + 2 J 


V4 


x bE^yy) 
3 n 


dx. 


\ vt lx) 

and for the //-parallel polarization we obtain 


Hi(x',f(x’)) - H(x'j(x’)) - exp(ikx' sin0 ( ) 

Gj w^['+ mr ]*"<**■> 


y=M 

y’=f{x') 


( 11 ) 


| dx ( 12 ) 


In Eqs. (II) and (12) the unknown quantities are bE/bn and H. Thus, these two integral 
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equations can be solved for either bE/bn or H on the surface of the grating. These quantities 
will be needed later in the analysis. 

The Green's function (i(x.v; x'.y') appearing in Eqs. (5) and (7) may be represented in 
various forms. The one form particularly suited to our needs is 


0(x,y,x',y’) = — 


exp {—ikx sindj) 
2ia 


m— — ° 


exp[— i27rm(x-- x , )/a]exp(+/k cosfl^ |y- 
k cos 8, 


-y'\) 


•m 


(13) 


Also, the diffracted field can be represented as a discrete angular spectrum of plane waves 
Ejix’.y') = ^ A m exp[i7c(x' sinO m + / cosi^)]. (14) 


In both Eqs. (13) and (14) the index (m) labels the diffraction order. Equation (14) de- 
scribes a plane wave whose wave normal makes an angle 0 m with the y axis and where the 
time dependence is given by the complex exponential exp(— /cut). 

It is important to understand that the expansion coefficients A m are constants (indepen- 
dent of x',y') only in the region away from the grooves of the grating. Inside the grooves 
they depend on both x and y’ as has been shown by Zaki and Kalhor. In the approximate 
theory of vector diffraction these coefficients were assumed to be constant everywhere. 

We will use the expansion of Eq. (14) for E d only in the region above the grating grooves. 
When Eqs. (13) and (14) are used in Eq. (5) we find the following expressions for the coeffi- 
cients A,„ 


X 



,A a£, (*>T) 

bn 


dx 


(15) 


for the E -parallel case. For the //-parallel case we obtain 


)/ mx - y) [-{^) s + 1;] 


X exp[-ik(x sin0,„ +y cos0 w )] 


y=f(.x) 


dx 


(16) 
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where B m is the expansion coefficient appearing in the discrete angular spectrum representa- 
tion of II j 

oo 

Ifj(x'y) = ^ B m exp[f&(x' sin0 m + y' cosO m )]. (17) 

m -— 00 

In both Eqs. (15) and (16) the total fields E and H appear in the integrands. These are the 
values that these fields assume on the surface of the grating profile. The angular spectrum 
expansions with constant coefficients given in Eqs. (14) and (17) cannot be used for them. 
These fields may be found by solving the integral equations given in Eqs. (1 1) and (12) as 
mentioned previously. 

In summary, the diffraction amplitudes A m and B m are obtained as follows. We first 
solve the integral equations given in Eqs. (1 1) and (12) to obtain the fields bE/bn and H on 
the surface of the grating. These quantities are then used in Eqs. (15) and (16) to determine 
the amplitude coefficients A m and B m 
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3. MAXIMUM DIFFRACTION EFFICIENCY USING 
RIGOROUS DIFFRACTION THEORY 

As reported in the 1971-1972 final report the problem of maximum diffraction efficiency 
was studied by using the computer programs based on the approximate theory. The amount 
of energy in a preassigned diffraction order was maximized by optimizing the initial solu- 
tion. It was also shown that the values of 0, and \/a required to maximize the solution were 
the same conditions required to produce a double Wood’s anomaly. We refer the reader to 
Fig. 5.3, p. 26, of that final report. We reproduce that curve here, in Fig. 3, for ease of 
reference. The plot shown in Fig. 3 is for the ^-parallel state of polarization. However, the 
figure is not complete because the calculation was not performed for the //-parallel polariza- 
tion. This is because the approximate theory was not able to produce reliable results for the 
//-parallel case. The calculation was later repeated using the Berkeley programs for both 
states of polarization. Below, in Fig. 4, we have reproduced the results of this calculation as 
they were first presented in the December 1972 monthly report. From Fig. 4 we see that 



Fig. 3. Normalized irradiance X/<z for the —1st order versus the ratio tya for the E- 
parallel state of polarization. The angle of incidence is —19.47122. 
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the calculation based on the rigorous theory corroborates the results obtained from the 
approximate theory for the ^-parallel case. In Fig. 4b we see that the efficiency curve for 
the //-parallel case is a maximum at the same value of X/a as was obtained for the /^-parallel 
polarization. It is also apparent that the irradiance in the —1st order quickly drops to lower 
values on either side of the blaze condition X/a = 0.6667 for the //-parallel case. Thus, we 
see that the grating is well blazed but only at a single wavelength. 

The occurrence of high blaze efficiencies simultaneously with the phenomenon of double 
Wood’s anomalies has also been reported by McPhedran. 3 In our studies this condition was 
obtained by optimizing the solution predicted by a geometrical model as reported in the 
monthly report dated 1 1 August 1972. At that time we reported on a graphical method for 
determining the grating parameters required to produce single and double Wood’s anomalies. 
That method is presented below and we further show that the conditions required to pro- 
duce a double Wood’s anomaly are also the conditions necessary for one of the diffracted 
orders to be diffracted back along the direction of the incident beam. In fact, at the condi- 
tions that produce a double Wood’s anomaly there is always one beam diffracted in the 
direction of the incident beam. 

From the grating equation 

sin0; + sin0 m = ~- 

where 0 ( - is the angle of incidence, Q m is the diffraction angle for the mth order, m is the 
diffraction order, X is the wavelength, and a is the grating constant, the ratio X/a is solved 
for and plotted as a function of the angle of incidence for a given order m. The value of 
sin0 m is set equal to ± 1 .0 because 9 m = ±90.0° when the mth order is evanescent. The curve 
so produced defines the values of X, a, and 0, required to produce a single Wood’s anomaly. 
The intersection of two curves defines the conditions required to produce a double Wood’s 
anomaly. Also, curves that lie above a given point on the graph indicate the real or homo- 
geneous orders present for a given set of parameters. Curves below a given point represent 
inhomogeneous or evanescent orders. A plot showing the conditions required to produce 
single and double Wood’s anomalies is shown in Fig. 5. 

If we now take the grating equation and set d m = 6,- and plot the ratio X/a as a function 
of 0,- for a fixed order m, we obtain Fig. 6. These curves then give the conditions necessary 
for the mth order to be diffracted back along the direction of the incident beam. Under 
these conditions the mth order is in autocollimation (Littrow mount). 

If we now superimpose Figs. 5 and 6 we obtain Fig. 7. This graph clearly shows that at 
the conditions required to produce a double Wood’s anomaly there is always one diffraction 
order that is autocollimated. It is also apparent that there are more than one set of condi- 
tions for which a given order can be used in autocollimation in the presence of a double 



C7> .£> 
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Wood’s anomaly. From these graphs the conditions required to obtain a high blaze effi- 
ciency for a grating operated in the autocollimation mode can be easily obtained. We are 
currently preparing a paper for the Journal of the Optical Society of America reporting 
these findings. 
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4. MAXIMUM DIFFRACTION EFFICIENCY - AN ANALYTIC SOLUTION 

The design of gratings that operate with a high blaze efficiency can also be persued ana- 
lytically. This approach requires that the equations derived from vector diffraction theory 
be solved for the grating profile fix) that will produce the desired blazed conditions. This 
procedure would be extremely difficult to do using the equations derived from the rigorous 
theory. Thus, the equations obtained from the approximate theory are used, and the results 
arc checked to insure that the fix) so obtained is consistent with the rigorous theory. In 
deriving the expression for the grating profile fix) we impose the constraint that all of the 
amplitude coefficients be exactly zero except one, the preselected blaze order denoted as 


the pth order. Thus, 

A m = 8 mp for all m. (18) 

The grating profile obtained from this study is given by 

fix) = ~[p(\/a)l2 cos9j]x. (19) 

The blaze angle is given by 

a = arc tan[— p{\/a)!2 cos0,] . (20) 

The angle of incidence is 

Oj = arc sin[p/(X/2a)] , (21) 

and the apex angle ^ is given by 

\p = 90.0° - a. (22) 


This grating profile and the related parameters are shown in Fig. 8. 



Fig. 8. Sketch of the profile f(x) for maximum blaze efficiency 
as predicted by theory. 
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The predicted grating profile shown in Fig. 8 was tested using the computer programs 
based on the approximate theory. The results indicated that the profile produced a grating 
that was perfectly blazed (100%) for both states of polarization in a preselected diffraction 
order when the grating was used in autocollimation. Furthermore, the grating was blazed 
over a wavelength interval of approximately 18.0 nm. The groove depths used in the above 
study varied from less than one wavelength to 3.3 wavelengths. However, the energy check 
revealed that even for the deepest grating the computer programs based on the approximate 
theory were operating perfectly. This result suggests that the approximate theory had 
allowed us to predict a grating profile for which the amplitude coefficients were not a func- 
tion of position within the grating grooves. This surprising result suggests that we (1) test 
the above profile with the Berkeley programs based on the exact theory and (2) investigate 
whether the above profile is also a solution of the exact equations when subjected to the 
same constraints. 

In trying to test the above profile on the Berkeley programs we found that these pro- 
grams would not run with the special profile and the constraint of autocollimation. How- 
ever, by varying the angle of incidence slightly, the Berkeley programs would run. The 
energy check revealed that the rigorous programs were not operating as well as might be 
expected. Under these conditions the rigorous programs indicated that 92% of the incident 
energy was contained in the preselected order (the —4th order was chosen) for the ^parallel 
state of polarization. It is believed that if the Berkeley programs had been able to operate 
under the ideal conditions, an energy check closer to unity would have been obtained to- 
gether with more energy appearing in the desired order. 

At this point we must mention that we have often encountered difficulty in getting the 
Berkeley programs to operate in a number of cases. Also, these programs and the approxi- 
mate programs are not able to analyze some of the gratings in which NASA-GSFC has ex- 
pressed interest. For these reasons we have embarked on a new method for obtaining a 
solution to the rigorous equations. In Section 5 of this report we describe these efforts in 
some detail. 

We now focus our attention on the investigation of whether the profile /(x) of Fig. 8 is 
also a solution of the exact equations. We start with the condition given in Eq. (18) and 
slightly modify it to read 

A m = 5 m p all m and y>0. (23) 

This modification sets a severe restriction on the behavior of the amplitude coefficients in 
the spaces above and within the grating grooves. It allows for no variation of A m as a func- 
tion of position within the grating grooves. 

If the conditions expressed in Eq. (23) are used, then the diffracted fields can be ex- 
pressed as an angular spectrum of plane waves with amplitude coefficients that are constant 
everywhere. The diffracted fields are therefore given by Eqs. (14) and ( 17). This means that 
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in the equations relating the amplitude coefficients A m and B m with the fields on the sur- 
face of the grating (Eqs. ( 1 5) and ( 1 6)) we may use the angular spectrum representation for 
the fields under the integral sign. Under these conditions the rigorous equations become 
much simpler, and we have shown (see Appendix A) that the special grating profile given by 
Eq. ( 19) is also predicted by the equations of the exact theory. We intend to publish these 
findings in the Journal of the Optical Society of America. 
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5. A NEW METHOD FOR THE SOLUTION OF THE 
RIGOROUS THEORY EQUATIONS 

One of the gratings that NASA-GSFC has expressed an interest in is the coarse echolle 
grating used in the vacuum ultraviolet. In general, these gratings tend to be very deep with 
respect to the wavelength of the radiant energy. Also, it is not uncommon for more than 
100 real diffraction orders to be present for a coarse echelle as opposed to only a few real 
orders as is usually the case for an echelette grating. As a result, it is not possible to analyze 
coarse echelle gratings by using existing computer programs. In addition, it is not possible to 
simply increase the number of diffraction orders that current programs can analyze because 
the programs soon become too large for the computer (CDC 6400). Thus, a new method of 
analysis must be formulated for echelle gratings. Below we outline a new approach for cal- 
culating the amplitude coefficients by using rigorous vector diffraction theory. For this 
method we hope to obtain a more convenient set of equations that can eve n be us ed when 
there are a large, number of orders present. In Appendices B, C, and D a mathematical de- 
scription of important parameters can be found. 

Before discussing the method of solution, we want to point out some of the difficulties to 
be expected in applying the rigorous vector diffraction equations to the study of echelle 
gratings. 

First, a computer can perform operations only on discrete or sampled values of a func- 
tion. The number of sampled values required to adequately represent a given function de- 
pends on how quickly the function varies: slowly varying functions require fewer sampled 
values than rapidly varying functions. Therefore, a function of two independent variables, 
k(x,x'), may require a very large two-dimensional matrix to store all of the sampled values 
required to represent k{x,x') if the function is rapidly varying. 

In the echelle problem for the £-parallel polarization, we wish to solve Eq. (11) for the 
unknown quantity d£(x,/(x))/dn. This quantity (matrix of sampled values) is then used in 
Eq. (15) to calculate the unknown amplitude coefficients for the diffracted fields. In order 
to solve Eq. (11) for the unknown quantity bE(x,f(x))/dn we must first calculate the ele- 
ments of the matrix representing the Green’s function G(x,/(x); x',/(x')) as given by Eq. 
(13). However, in Eq (13) the index m is the diffraction order and appears in one of the 
exponentials. For large values of m (up to 100) this exponential will be a rapidly oscillating 
function and, therefore, the Green’s function will require a very large matrix to contain all 
of the necessary sampled values. A simple calculation shows that for 1 00 real orders the size 
of this matrix alone exceeds the available core storage of many large computers. 

One way of overcoming the above storage requirement is to express the Green s function 
in terms of functions that are more slowly varying than the plane wave representation Eq. 
(13). Another possible representation for the Green’s function is in terms of Hankel func- 
tions. This representation converges rapidly for gratings with a period much greater than the 
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wavelength of the radiant energy. This condition is easily met with coarse echelle gratings 
but not necessarily with echelette and holographic gratings. Therefore, the Hankel function 
representation is of limited use in our work. 

It is not unusual for functions that are rapidly varying in the spatial domain to be much 
smoother in the spatial frequency domain. For such functions the analysis can be performed 
in the frequency domain by using smaller matrices. In the work presented below we have 
chosen to perform the calculations in the frequency domain with the aim of reducing the 
sizes of the required matrices. This is done by expressing all of the quantities in terms of 
their complex Fourier series expansions. This method has an additional benefit due to the 
uniqueness of the Fourier series. Namely, that the plane wave and the Hankel function 
representations of the Green’s function (as well as any other representations) must have the 
same Fourier representation. This means that the Fourier method can be used for the analy- 
sis of echelcttes as well as coarse echelle gratings if the matrices are of reasonable dimensions 
in the case of the latter. 

To illustrate the proposed method of analysis we note that Eq. (1 1) is of the form 

h{x’) =( k(x,x')g(x)dx, (24) 

o 


where 


/»(*') = Efix’JKx’)) 

k(x,x') = / exp(+/kx' sin0,)G(x/(x);x\/(jc')) 


(25) 


and 


Ii(xJ'(x)) = k\(xj(x)) + E d (xJ\x)). 

We then expand each of these functions in terms of their Fourier series representations 

o© 

/?(x r ) s h k expi+ilirkx' /a) 


k(x.x') 


g(x) = 


ft 




kj,, exp(-i2nnx/u) <zxp(+i2njx'/a) 


— — <» — oo 

gi exp(+i2nlx/a). 


(26) 
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The complex expansion coefficients are then given by 
1 f a 

h k = — I h(.x')cxp(-i2Trkx'/a)dx' 
aJ Q 

i f-a ft* 

kj n ~ ~2 I I k(x,x') exp(+i2imx/a) exp(-i2irjx'/a)dxdx' (27) 

a J 0 } 0 

1 ( a 

g t = - I g(x) exp(-i2irlx/a)dx. 
a 0 

Substitution of Eq. (26) into Eq. (24) yields 

no « tro 

y h k exp(+i2nkx'/a) =f ^ ^ k jn exp(-i2irnx/a) 

0 n =—°° /=- 00 


DO 

X exp(+/27r/x7a) ^ g; exp(+/2irbc/tf)fl?x. 

/=— 00 


(28) 


After some manipulation Eq. (28) becomes 

<*> 00 / °° \ 

V (A fc )exp(+i27r*x7a) = ^ [a ^ exp(+/27r/x7a). (29) 

h-~\ «=— / 

For Eq. (29) to be valid for all terms in the expansion we must have 



Equation (30) is a matrix equation where h and g are column matrices and k is a two- 
dimensional matrix. The size of the matrices is determined by the values of; and n beyond 
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which the contributions of the series for h and k are less than the specified error. That is, n 
and j cannot assume the infinite limits expressed above but must be limited to some maxi- 
mum value. The maximum values are chosen such that the errors due to neglecting terms 
beyond these maximum values are less than a specified amount. 

We wish to determine the unknown expansion coefficient g n in Eq. (30). This can be 
done because we can calculate the coefficients hj and k: n because they involve known quan- 
tities. In Appendix B we have derived the equations from which the k ]lt coefficients are to 
be calculated. In Appendix C we derive the equations for the coefficients hs. 

Note that / and n can assume negative as well as positive values. If we denote the mini- 
mum values as JMIN and NMIN and the maximum values as JMAX and NMAX then the 
number of values that j and n will assume are given by J and N, respectively, where, 

./ = JMAX - JMIN + 1 

N = NMAX - NMIN +1. (31) 

In general, we will have J = N . This means that we can regard Eq. (30) as either a matrix 
equation or as N simultaneous equations. Solution of the matrix equation will require that 
the inverse matrix kjj, be calculated. In any case, once the coefficients g„ are solved for the 
function, g(x) is determined. Note that g(x), as defined in Eq. (25), also appears as the 
unknown in the equation for the amplitude coefficients Eq. ( 1 5), namely, 

— )f expl-ik(x sin6 m +f(x) cosd m ]g(x)dx. (32) 

m / J o 

We now expand the complex exponential in terms of its Fourier series representation 




1 


2ak cos 9 


e m (x) = &xp[-ik(x sind m + f(x) cosd m )] 



e m n exp(— Hirnx/a). 


(33) 


Thus, 


^ m \2 ak eos0 


\ <1 j” * °° 

l ~ Q-)f J e m .n exp(-i2irnx/a) ^ g, exp(+i2nlx/a)dx 

“ m ' 0 /!=—■» /=— oo 


I 


2 ak cosd 


oo oo 

S 2 


JTl f}~ - oo /=— a 


nSl a ^l,n 
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Equation (34) is also a matrix equation where A andg are column matrices and € m n is a 
two-dimensional matrix. The coefficients e m „ are to be calculated from 

1 f a 

e mn ~ ~ e m (x)exp{+i2irnx/a)dx. (35) 

aJ 0 

Explicit expressions for the e m n coefficients are derived in Appendix D. Once these coeffi- 
cients have been calculated we can then calculate the unknown amplitude coefficients A m 
using simple matrix multiplication. 

We note that the method of analysis is the same for the /^-parallel state of polarization 
except that g(x) would now be proportional to the magnetic field and k(x.x') would repre- 
sent the normal derivative of the Green’s function. 
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APPENDIX A 

CONSISTENCY OF SPECIAL BLAZED GRATING PROFILE WITH 
EQUATIONS DERIVED FROM THE RIGOROUS THEORY 


The amplitude coefficients for each of the diffracted orders are given by 


A,n ' 2ak 


c cos 0 m ^ Iq 


exp [ —ik(x sin0 w + f(x) cos 6 )] 


* BEjxJjx)) 

Bn 


dx. 


We assume a plane wave expansion for the diffracted fields of the form 


oo 

^ A m exp[/A'(x sin0 OT + y cos Q m ) | . 




We <i Iso note that 


9 

dn 


WuA ± . A] f, + 

dx ) dx By] |_ l dx ) 




(AI) 


(A 2) 


(A3) 


If we use Eq. (A3) to operate on the total electric field and we substitute into Eq. (Al), we 
obtain after considerable manipulation 


f a ( df(x ) . „ \ 

\ rfx sin( ^j + cosB A exp[— ikf(x)(cosd n + cos0,)] exp(—i2nnx/a)dx 


S A rn | jf {~^ sin ^OT - cos #mJ exp [ —ifcf(x )(cos0„ - cosf) m )] 


X exp[— /2 tt(h - m)x/a]dx + 2a cos0 m 5 m „ | . 

In deriving the equation of the grating profile we imposed the constraints 

A m ~ *>mp’ for all m 
0 m = Or 


(A4) 


(A5) 
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These constraints are then used to simplify Eq. (A4). The left-hand side (LHS) of the equa- 
tion then becomes 


LHS = - sin#.- f exp(-i'2kf(x) cos#; - i2irpx/a)dx 

' J 0 dx 

c a 

- cos#,- J exp(-iW(x) cos#,- - i2itpxfa)dx. 

"0 

The right-hand side (RHS) reduces to 
RHS = f(a) sin#,- +« cos#,-. 

When we match coefficients for sin#, in Eqs. (A6) and (A7) we obtain 


(A6) 


(A7) 


f(a) = - f exp(— i'2fc/(x) cos#,- - i2itpx/a)dx. 

J o ax 


(A8) 


The grating profile /(x) is given by 


Ax) = 


/ pW V . 

\2 cos#,- / 


(A9) 


Substitution of Eq. (A9) into Eq. (A8) for the value x = a yields the identity f(a) - /(a). 
Thus we see that Eq. (A9) is a solution of the rigorous equations under the constraints 
expressed by Eqs. (A2) and (A5). 
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APPENDIX B 

EXPANSION COEFFICIENTS k jn 
In Section 5 of the text we defined 

k(x.x') = i expl+ikx' sin©, ] G( jc/(x); x',/(x')) (Bl) 

where Gix/ix); x'J(x')) is the Green’s function given by 


G(xJ(x);x',f(x’)) = - 


exp(+iftx sin#,-) 

2ia 


X 





m 


exp[-/2ffw(jr - x')/a] exp (+ift cosg ffl j/~(x) -/(x')l] 

ft COS0 OT 


(B2) 


Recall that we expressed kix,x') in terms of its Fourier series expansion 


ft(x.x') = 2 2 kj n exp(— i2imx/a) exp(+i27r/V/a). (B3) 

n-—°° /=—«• 

Therefore, the expansion coefficients are given by 


In order for these coefficients to be calculated we insert Eq. (B2) into Eq. (Bl) and the 
resulting expression is used in Eq. (B4), 

In order for the calculation to be completed we must know the equation of the grating 
profile fix). 

For this derivation we assume fix) to be the profile shown in Fig. Bl. The actual grating 
profile is given by 


ft(*„x') exp(+/27rnx/a) expi-ilirjx' la)dxdx' . (B4) 

0 



Ax. 


Ax± 
a — x x 


(a - x). 


fix) = 


0 < x < x x 


(B5) 
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Fig. B1 . Grating profile f{x). 


where A is the slope of /(x) in Region I ( A = tana). After much calculation we find the 
expression on page 26 for the Fourier expansion coefficients k } - n for the profile shown in 

Fig. Bl. 
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K j,„ *1 (V 4 ’' 2 )' 


P 0™0 

12 


m=-°° 


a- jP 


to - 

12 


; + i2TrCn-j)x 1 /a _ e + i2 ir Cl/pffl ~ j/a)x-, 

Ta-nPj“ J )/pW 


(m) 


11 \f e - i2 ^(l/ p ii -n/a)x 1 _ x e + i2Tr(n-j) Xl /a 


a-jP 


jo" TO] 


11 


- IkBaCosO 
e m 


(a-nP 1 ^Vp[ I J 0 


a pCm) 

P 12 


n(™) 
22 

V-^IF 

(m) 


(n-j ) 

+ i2TrCa/P^2' ) - j) ,-i2TTY 1 x 1 


1 ) 


Va-nP 


TO 

12 


Wi2TrU/P 2 ^TO/a) X;L ^ e -i2Tr(l/P 


1 ( 2 ) -n/a)x 1 . x 
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where. 




m 

a 



m 

a 


21 + 


a 


ACos6 w 

m 

A 


) 


ACos9 m 

m 

X 


) 


BCosO 

m 

X 

BCosQ 

m 

X 


) 

) 


A = Tan a, B = -Tan 3 
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APPENDIX C 

EXPANSION COEFFICIENTS h } 

Using Eq. (27) of the text for the expansion coefficients hj we have 

1 f a 

hj - -J h(x') exp\~i2irjx'fa]dx' 
a J o 

where h(x') is the incident electric field 


(Cl) 


//(*') = exp [—/Ay sinO,] expl-tfr/t*') cos0,]. (C 2 ) 

We use the same grating profile /(x') defined in Eq. (B5) of Appendix B. In this way we find 
the following expression for the Fourier expansion coefficients 

h i " (i){a [expC-^TroorO- U + (^)exp(-/A:&t cosfl,) 


where 



(C3) 


« = (sin0,/A + A cosOj/X + j/a) 
fi = (sin0,/X - B cos Oj/X + j/a). 
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APPENDIX D 

EXPANSION COEFFICIENTS e m n 

Using Eq. (35) of the text for the expansion coefficients e m n we have 

i c a 

e m,n ~ -] e m(x) zxp(+i2Tmx/a]dx (Dl) 

aJ o 

where 

^,,(*) s exp[-/*(xsin0 m +/(*)cos0,„)]. (D2) 

We use the same grating profile /(x) defined in Eq. (B5) of Appendix B. In this way we 
obtain the following expression for the Fourier expansion coefficients 

e m,n = ( 2 ^) {(^) [exp(-/2jrax, ) - 1 J + exp(-i2nBa cos 8 m 

X [exp(~i27rj3<z) - exp(-f27rj3x, )] | (D3) 

where 

a: = (sin0 m /X + /l cos 0 m j\ ~ n/a) 

P = (sin0„, /X - B cos0 m /X - n/a). 
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